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DESIGN  AND  MAXIMUM  ERROR  ESTIMATION  FOR 
SMALL  ERROR  LOW  PASS  FILTERS 

by 

W.D.  Hibler  III 


Introduction 

There  exist  many  situations  in  which  it  is  desirable  to  filter  out  certain  frequency  components 
from  geophysical  time  series  with  a  sharp  frequency  cutoff.  When  using  a  finite  nonreenrsive  digital 
filter  for  this  purpose,  some  smoothing  of  the  frequency  cutoff  must  be  employed  to  reduce  the  magni¬ 
tude  of  spurious  frequency  ripples  in  both  the  stop  band  and  the  pass  band.  One  type  of  smoothing 
for  this  purpose  was  developed  by  Martin  (1957).  In  later  work  numerical  search  techniques  for 
finding  the  optimal  type  of  smoothing  in  terms  of  minimizing  the  ripple  error  for  a  given  transition 
bandwidth  were  discussed  by  Gold  and  Jordan  (1969).  However,  in  both  these  techniques  the  design 
procedure  depends  upon  estimating  errors  numerically  for  individual  filters.  For  convenience  it 
would  be  quite  useful  to  have  a  design  procedure  allowing  the  accurate  estimation  of  ripple  error  for 
any  low  pass  filter  designed  by  using  the  procedure.  This  is  especially  important  for  filters  con¬ 
sisting  of  large  numbers  of  weights  because  the  error  calculation  is  time  consuming  and  cumbersome 
in  such  cases. 

Such  a  formulation  and  an  error  estimation  are  carried  out  here.  The  smoothing  is  carried  out 
by  using  standard  spectral  windows,  developed  for  purposes  of  spectral  analysis  (Blackman  and 
Tukey  1958),  Specific  examples  illustrated  later  show  that  such  a  smoothing  procedure  yields  low 
pass  filters  with  smaller  errors  (by  a  factor  of  five  in  one  case)  than  those  of  comparable  filters 
designed  by  the  Martin  (1957)  technique.  As  a  result  of  the  error  estimation,  the  straightforward 
working  equations  may  be  used  to  quickly  calculate  convolution  weights  for  symmetric,  unity  gain 
low  pass  filters  with  a  guaranteed  ripple  amplitude  outside  the  transition  band  of  less  than  0.9% 
or  0.05%  depending  on  how  wide  the  transition  band  is  made.  Moreover,  in  the  working  equations 
presented  later  the  width  of  the  transition  band  scales  inversely  with  N  (the  number  of  convolution 
weights),  so  that  the  frequency  cutoff  for  a  given  error  may  be  made  arbitrarily  sharp  by  increasing 
N.  Once  the  convolution  weights  are  calculated,  the  convolution  may  be  performed  by  conventional 
techniques  or  by  the  aperiodic  fast  Fourier  transform  (FFT)  procedure  (Stockham  1966),  a  much  more 
rapid  technique.  (The  aperiodic  FFT  procedure  should  not  be  confused  with  simply  doing  a  FFT, 
removing  certain  frequencies  and  then  transforming  back  to  real  space,  because  such  an  operation 
effects  a  periodic  convolution  with  spurious  end  point  effects).  For  completeness  we  have  included 
the  basic  equations  for  the  FFT  aperiodic  convolution  procedure  in  Appendix  A. 

In  the  next  section  we  shall  formulate  the  design  procedure  in  a  manner  that  allows  the  numerical 
estimate  of  maximum  error.  Maximum  errors  will  then  be  calculated  for  the  three  most  standard  spec¬ 
tral  windows.  A  specific  example  and  working  equations  will  also  be  presented  for  those  interested 
in  using  low  pass  filters  without  going  through  the  complete  detail  of  the  design  procedure. 
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Filter  design  and  error  estimation 

A  convenient  way  of  viewing  the  frequency  response  of  a  given  finite  filter  is  to  consider  the 
frequency  spectrum  of  the  filter  an  “estimate”  of  the  true  spectrum  of  the  infinite  length  filter. 
Viewed  in  this  manner,  spurious  frequency  ripples  result  from  the  convolution  of  the  Fourier  trans¬ 
form  of  the  data  window  with  the  true  spectrum.  A  way  to  improve  the  situation  is  to  use  a  uon- 
rectangular  data  window  whicli  produces  smaller  amplitude  ripples,  but  which  also  smooths  out  the 
frequency  cutoff. 

Let  us  consider  a  symmetric  digital  filter  with  weights  C(n ): 
n  0,  +  1 . +  Af  -  1 

n  =  *  N. 

The  frequency  response  of  this  filter  is  (Holloway  1958): 

AT  1 

HU)  2  1  H(n)  cos  (2n(n\  t)  i  H( 0)  i  cos  (2n(NN  t)H(N)  (1) 

u  t 


C(n) 


H(  a) 


H(n ) 


where  At  is  the  data  interval.  We  shall  denote  functions  in  frequency  space  by  a  tilda  variable 
witli  I  ranging  from  0  to  0.5  cycles/data  interval.  If  R{t)  is  known  at  the  discrete  frequencies  n/(2T) 
where  T  N At),  then  eq  1  may  be  inverted 

.  1  nkn  G,,  /o<m  W(0)  ,  .  H(N/2T) 

H(n)  -  >  cos  — -  H(k/2T)  h  -L-  i  cos  ( m )  ■  ----  .  (2) 

N  N  2N  2 N 

k  1 

Introducing  the  infinite  Dirac  comb 


oo 

V(r,  At)  1  «(t  -  qA t) 

q=— oj 


(3) 


(where  fi(l)  is  the  Dirac  delta  function),  eq  1  may  be  rewritten  in  the  form 


HU) 


dl'  cos  (2ntv)W(i')K](r\  At)  /  H^V  -  f  j)V(t 2T)dil 


(4) 


where 


WU) 


n^t)  = 


o  let 

->  T 

1/2  \V\ 

T 

i  in 

v  T 

/7(n) 

n 

o, . 

H(n)/ 2 

n 

*  N 

0 

HI  " 

T. 

(The  values  of  //jft)  at  other  values  of  l  are  not  important  because  of  the  multiplication  by  V(f  ;  At) 
in  eq  4.  To  verify  that  eq  4  is  iu  fact  equal  to  eq  1,  it  is  convenient  to  note  that  the  convolution  of 
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ite  Dirac  comb  is  equal  to  £  Hj(t'  -  2gT)  so  that  foi  \t  |  ^  NAt  this  sum 


q=-° 


H  j(t)  vyith  the  infinite 
over  q  reduces  to  H  j(t' ).) 

Since  all  terms  in  the  integral  in  eq  4  are  even  functions  of  t\  we  may  replace  cos  (2n(P)  by 
ei2nft'  Using  the  identity 


;  A  t)dt'  -  —  V(f;  1/At) 


re|M'  v(f 

— no 

and  the  convolution  theorem,  eq  4  becomes 


(5) 


H(0  =  l  I  W(f  -  tt/2 T  -  q/At)H(n/2T) 

n-  q=~oo 


(6) 


where 


m  .  -■*"  .  «„(/). 


(7) 


Eouation  6  illustrates  how  the  aliased  Fourier  transform  of  the  spectral  wtndow  “ 

the  discrete  frequency  response  to  give  the  complete  frequency  response  We 

change  the  form  of  W(f)  by  using  some  nonrectangular  spectral  window, call  it  *  V W. 

of  the  interval  t-T.  T).  If  we  substitute  such  a  data  window  m  eq  4,  then  eq  6  is  the  same  . 

W(()  being  replaced  by  W^t)  and  H(0  being  a  different  frequency  response,  call  it  H 

We  would  like  to  use  a  spectral  window  W^t),  having  as  small  side  lobes  as  possible  while 
causing  'minimum  smoothing  of  the  spectral  content.  A  general  class  of  windows  developed  for 
power  spectral  calculations  with  these  considerations  in  mind  is  of  the  torm 


ai0  +  2£aijcosT- 

it!  <  T 

j=i 

0.5[aiO  +  2  1  (-1)4  1 

Itl  = T 

j=l 

0 

Itl  > T ■ 

(S) 


The  Fourier  transform  of  any 


window  of  this  type  is  given  by  (Blackman  and  Tukey  1958,  p.  99) 


QXO  =  ai0Qi0(/)  +  S  a^lQoCf  +  j/2T)  +  Q0(f  -  j/2T)] 
j=1 

where  Qn(f)  is  as  defined  in  eq  7.  In  particular,  we  will  consider  only  cases  where  a, ,  =  0  for 
T>  2  In  this  case,  since  Q0(m/2T)  =  0  for  an  integer  m  with  magnitude  >  0,  we  have 


I  Qj(n/2T  -  Nq/T)  =  Q j(n/2T) 

q=r-00 


(10) 
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for  (n i  •-  N  -  3.  Consequently,  eq  8  [with  WU)  replaced  by  W^f)!  takes  on  the  simple  form  at 
discrete  frequency  values  n/2T. 

in  2  '  ' 


it  H(n/2T)  0  for  |n|  N  3.  Equation  11  can  thus  be  used  to  determine  smoothed  discrete 

frequency  response  values  H  ,<n/2T)  {riven  ideal  frequency  response  values  H(n'2T)  and  some 
spectral  window  coefficients  ajn).  Filter  weights  may  then  lie  calculated  from  eq  2. 


We  will  consider  here  die  three  most  standard  spectral  windows  (Blackman  and  Tukey  1958) 
which  have  die  listed  nonzero  coefficients  aj(||: 

1.  Hanning  a,  0 

0.5,  a,  ,  a  j  ( | 

0.25 

2.  Hamming  aL,  0 

0.54.  ag<  j  u 2  ,  ] 

0.23 

3.  Blackman  a;j  Q 

0.42,  a, ?,  j  a3>  ( , 

-  0.25,  a3  r  a3  2  0.04  . 

Using  eq  11  the  design  procedure  is  to  choose,  initially  IUn/2T)  equal  to  1  up  to  [n|  p,  and 
zero  for  jn|  p.  We  will  call  p  the  design  cutoff.  Smoothed  frequency  weights  W,(n  '2T)  are 
determined  by  eq  11  and  the  smoothed  frequency  weights  are  converted  to  digital  filter  weights  by 
eq2.  Clearly  the  transition  hand  will  extend  from  (p  -  1) '2T  to  (p  ,  2)/2T  for  the  Hamming  and 
Hanning  spectral  windows  and  from  (p  -  2),  2 T  to  (p  .  3V2T  for  the  Blackman  window.  Errors  are 
any  deviation  from  1  in  the  pass  band  and  deviations  from  0  in  die  stop  band. 

The  basic  equation  for  estimating  the  maximum  error  is  eq  8.  An  equivalent  form  of  this  equation 
obtained  by  rearranging  the  double  sum  is  (frequencies  are  measured  in  units  of  1/2T): 


«!</)  1  -  [)  (12) 


where //A<f)  1  H((  2qN)  is  now  periodic.  Keeping  the  periodicity  of  WA(/)  in  mind,  we  see 

(j  >•  A 

that  if  we  search  numerically  for  the  maximum  error  we  need  consider  only  low  pass  filters  witli 
design  cutolt  values  p  ^  N  2.  This  follows  because  a  filter  witli  a  higher  design  cutoff  creates  the 
same  sawtooth  function  (witli  a  shift  in  the  origin)  for  WAU)  as  for  1  WA'(/),  where  //A'(/)  is  an 

appropriate  set  of  design  weights  witli  design  cutoff  p  less  than  N-' 2.  But  the  maximum  error  for 
1  «A*</)  is  ,,1<>  SiII,le  as  fo1  wa,(/)  Because  convolving  WU)  with  1  yields  1  identically.  (As 

before,  by  maximum  error  we  mean  the  greatest  deviation  from  1  in  die  pass  band  or  from  zero  in  die 
stop  band,  whichever  is  larger.  Thus,  changing  the  origin  of  a  given  sawtooth  function  does  not 
change  the  maximum  error  defined  in  this  way.) 

Now  for  design  cutoffs  less  than  or  equal  to  N/2  we  may  write  eq  12  in  die  form 


inn 


o 

1  W  (k)H.(k  -  f)  +  E 
k  -N(, 


(13) 


when1 
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£  "  ,  J  .  A <k  -  I)  *  1  WlkW.tk  . 


k=N0+l 


k=-N0-l 


/) 


(14) 


the  maximum  error  introduced  by  negS^  of  thefwm1^)6 

sin  nk  where  g(k)  is  a  monotonically  decreasing  function  past  a  certain  value  of  k  For  example*  for 
the  Hamming  spectral  window  g(k)  is  given  by  l/n{0Mk2  -  0.54 /k(k2  -  1)1  T 


No 

1 

k=-N 


W(k)H(K  -  /)  =  |y(f  -  p)  +  w(l  -  P  +  1)  +  ...  +  W([  +  p) 


(15) 


(where  p  is  the  design  cutoff  value)  for  value*  nf  n  _  1  s>  u 

calculation  was  performed  numerically  by  evaluating  ea  l5  at’th  °f  1  “P  t0  N°'  SuCh  a 

1/20  2/20  N  fnr  value*  ^  o  evaluatmS  e(J  15  at  the  discrete  frequency  points  t  =  0, 

absolute  value  of  ?//(/)  Ilf  r  ■  J  . N°'  The  maximum  error  was  determined  by  taking  the 

r"  ”;5,  !  '’  ,T  11  fW  '  “  the  pass  tand  al»  «»  absolute  value  of  HU)  fo,  ,  in  the  stop 

d.  S  nee  HU)  typically  goes  through  one  oscillation  for  A  f  =  1  [because  W([)  is  of  the  form  Mr) 

7c  2L‘ »  W  “"V4  ta  —  <—•««.  -  XmZJZSZS? 

of  change  in  HU).  Consequently,  errors  due  to  slightly  nUnning  the  maaimum  side  lobe  .ere 

1  T°  “"al““e  ,he  er™  'or  (liters  with  numbers  of  weights  2 N  +  1  less  than  (or  ential  to)  4N 
o’  1/20  “n  (untt^CtoVi/aT("  e”1Ual‘ng  lh0  ^OdOouoy  response  at  values  of  the  frequency  (  ° 
of  the  **  Vata  -  '  — ***  <•»*  -tao 

utrto  .ss-r was  be=r 

Uhd  0,  cutoff  values  boT^ll^  «* 

|  IB“K  ‘ ,he  nlter  orror  tor  N  greater  than  a  given  cutoff  value  is  illustrated  bv  a  straieht 

sf bi,i??T  •  £  rrsxsLS  ssr  z.. 

a  *&•'  Errors  using 

erras^My^temHhan8M4%^ "* 

Waiting  aquations 

~  rr:  “°w  - — 

=es  below  (N ,  /2Af)  cycles/data  K 
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Figure  1.  Maximum  error  versus  number  of  filter  weights  for  low  pass  digital 
filters  designed  using  standard  spectral  windows.  The  dotted  lines  denote  max¬ 
imum  error  calculated  accordir/g  to  eq  15  and  represent  an  upper  limit  that  is 
never  exceeded.  As  can  be  seen,  the  upper  limit  becomes  more  refined  for  large 
N.  For  a  given  N,  there  are  2N  +  1  digital  filter  weights. 


HU) 


1,  i  <  N  j 
0.77,  i  „  N j  4  1 
0.23,  /  =  Wj  +  2 
0,  j  s  Af,  ,2 


(16) 


for  the  Hamming  window,  and  frequency  weights: 


HU) 


0.91 

0.7 

0.2! 

0,0' 


<"i 

=  N\  4  J 
=  JV ,  f  2 
-  N,  +  3 
=  JVj  ^  4 
'  JVj  i  4 


(17) 


for  the  Blackman  window.  The  transition  band  frequencies  extend  (units  of  cycles/data  interval) 
from  N  j /2N  to  {N  j  +  3)/21V  for  the  Hamming  filter  and  from  Af  j/2Af  to  (Af  t  -i  5)/2A!  for  the  Blackman 
filter.  The  filter  weights  C(n)  are  obtained  from  eq  18 


1 

N 


At  —  1 


H(  0) 


//(Af) 


C'(n) 


2N 


2  N 


cos  irn 


(18) 
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Figure  2.  Amplitude  response  of  selected  nonrecursive  digital  filters  consisting  of  61  symmetric 
weights.  The  error  and  transition  bandwidth  (A  I)  in  cycles/data  interval  lor  the  various  filters  are: 
1.  (Hamming)  0.42%.  A /  =  0.050;  2.  (Martin)  0.62%,  At  -  0.054;  3.  (Blackman)  0.03%,  Al 
0.083;  ai m  4.  (Marlin)  0.17%,  Al  =  0.10.  For  the  Martin  filters  b  and  d,  the  rQ  parameter  is  0.234 
and  the  H  parameters  are  respectively  0.27,  0.051. 


where 

C(n)  -  C'(n)  for  n  =  0,  il,  +2 . +(N  -  1) 

and 

C(n)  -  C'(n)/2  for  n  =  iN. 

Note  that  for  numerical  computation  the  sum  over  i  in  eq  18  need  only  extend  up  to  the  last  nonzero 
value  of  H(i).  The  resulting  weights  C(n)  can  then  be  used  to  filter  some  time  series  rjii)  by  the 
convolution 

rjAi )  =  I  C(j)  t p  -  j)  (19) 

j*-N 

with  7]{(i)  being  the  filtered  result.  The  filtering  procedure  necessitates  the  loss  of  2 N  points.  As 
long  as  the  transition  band  is  contained  between  zero  and  the  Nyquist  frequency  (1/2  cycles/data 
interval),  the  maximum  error  is  less  than  0.9%  for  the  Hamming  filter  and  less  than  0.05%  for  the 
Blackman  filter  (Fig.  1). 

In  particular,  taking  N  =  30  and  (V  t  =  14,  Figure  2  illustrates  the  frequency  response  of  the 
filters  designed  using  eq  16-18.  We  have  also  illustrated  typical  filters  generated  according  to  the 


proceduredescribed  by  Martin  (1957).  The  working  equations  for  the  Martin  filters  are  also  given 
by  Davis  (1971),  a  somewhat  more  accessible  reference.  For  the  Martin  filters  there  are  two 
parameters  rc  and  H,  where  rc  denotes  the  100%  cutoff  frequency  and  H  determines  the  amount  of 
smearing  of  the  frequency  cutoff.  (Values  of  rc  and  H  are  given  with  the  figure  caption.)  Doth  the 
Hamming  and  Blackman  windows  generate  filters  with  less  side  lobe  error  than  the  Martin  filters 
witn  comparabie  transition  bands.  The  Blackman  window  has  decidedly  less  error  than  the  compara- 
ble  Martin  filter  (smaller  by  a  factor  of  5).  Moreover,  the  transition  band  is  considerably  smaller  in 
e  ackman  case,  with  a  width  equal  to  5/61  cycles/data  interval  as  opposed  to  about  6/61  cycles/ 
data  interval  for  the  Martin  filter.  A  Martin  filter  with  a  smaller  transition  bandwidth  would  be 
expected  to  have  an  even  greater  error. 


Application  example 


To  illustrate  the  application  of  a  low  pass  filter  to  digitized  data  we  have  applied  a  low  pass 
filter  to  laser  profilometer  data  taken  from  an  aircraft  flying  over  the  arctic  pack  ice.  The  results 
are  illustrated  in  Figure  3.  The  illustrated  data  consist  of  1800  digitized  points,  at  1-yd  intervals 
representing  a  straight-line  profile  of  the  upper  surface  of  the  ice  pack.  The  original  record  length 
was  2720  points  as  the  filtering  process  necessitated  the  loss  of  460  end  points  at  each  end  of  the 
record.  The  low  frequency  trends  in  the  data  are  due  to  the  aircraft’s  change  in  altitude  (about 
20  ft)  over  the  1800-yard  record.  The  smooth  curve  represents  the  low-pass  filtered  result. 


The  low  pass  filter  was  designed  using  the  Hamming  window  and  consisted  of  921  weights  with 
he  transition  band  extending  from  3/460  to  6/460  cycles/yd.  By  reference  to  Figure  1  the  error  is 
ess  than  0.6/„.  The  filter  weights  were  convolved  with  the  record  using  the  aperiodic  fast  Fourier 
transform  procedure  with  a  Fourier  series  length  of  2048  points.  In  this  particular  case,  the  FFT 
aperiodic  convolution  took  about  45  sec  (requiring  5  applications  of  the  FFT  algorithm  for  2048 
points),  which  was  about  five  times  faster  than  the  convolution  time  using  a  conventional  program. 


FvT'erl  EXan>P!e  °l  l°W  PaSS  'iltering-  The  i,,us£rated  da‘a  insist  of  1800  points  at  1-yd  inter- 

IZT  CT  iS  lhe  l0W  pass  tiIteKd  result  USing  3  Hamming pass  niter  consisting 
of  921  weights  with  a  transition  band  extending  from  3/920  to  6/920  cycles/ yd. 
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We  remark  that  the  low  pass  filter  gives  a  reasonable  estimate  of  the  aircraft  motion  (especially 
in  this  example),  but  not  in  general  ideal  because  there  is  an  overlap  between  the  surface  roughness 
spectrum  and  the  aircraft  motion  spectrum,  [A  more  detailed  filtering  process  which  bypasses  the 
spectral  overlap  problem  is  described  by  Hibler  (in  prep).]  Figure  3  is  a  good  example, 
however,  of  where  small  filter  errors  in  the  pass  band  are  important  because  of  the  much  larger 
amplitude  of  the  low  frequency  components  compared  with  high  frequency  components.  When,  tor 
example,  the  low  pass  filtered  result  is  subtracted  from  the  initial  profile,  large  errors  in  the  pass 
band  leave  residual  low  frequencies  with  amplitudes  commensurate  with  the  surface  roughness.  We 
note  in  passing  that  it  is  clear  that  rather  than  low  pass  filter  the  data  and  then  subtract  the 
filtered  points  the  same  end  result  could  be  obtained  by  applying  the  higli  pass  filter  consisting  of 
weights  Ct(n),  where  Cl(n)  -  S„t0  -  C(n)  with  C(n)  being  the  above  low  pass  filter  weights. 


Conclusion 

In  conclusion,  the  working  equations  described  in  eq  16-18  give  a  rapid  method  of  determining 
low  pass  filters  with  sharp  frequency  cutoffs  and  small  spurious  side  lobe  errors.  Moreover,  since 
the  maximum  filter  error  is  known,  individual  filter  errors  do  not  have  to  be  checked  -  a  laborious 
procedure  for  large  N.  In  practice,  the  procedure  generates  filters  with  small  errors,  especially  for 
the  Blackman  window,  which  are  smaller  than  those  of  comparable  filters  generated  by  existing 
techniques.  It.  is  probably  true  that  for  a  particular  cutoff  and  convolution  length  a  more  accurate 
nonrecursive  filter  could  be  designed,  for  example,  by  using  a  slightly  different  spectral  window. 
However,  in  most  applications  eq  16-18  give  adequate  filters;  and  these  equations  are  convenient  to 
use  because  a  maximum  filter  error  estimate  is  known. 
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appendix  a.  fft  aperiodic  convolution  technique 


Given  a  set  of  filter  weights  C(n), 
the  convolution 


we  would  like  to  filter  the  space  series  77(1)  by  performing 


Vf(i) 


N 

1  C(j)r,(i 
j--N 


;) 


(Al) 


wit  >7f(i)  being  the  filtered  result.  The  filtering  procedure  necessitates  the  loss  of  2N  end  points. 

or  a  filter  consisting  of  a  large  number  of  weights,  a  conventional  convolution  program  takes  " 
pro  lbitive  amounts  of  time.  Consequently,  it  is  necessary  to  use  the  aperiodic  fast  Fourier  trans¬ 
form  high-speed  convolution  technique  (Stockham  1966).  The  aperiodic  FFT  convolution  procedure 
should  not  be  confused  with  simply  carrying  out  a  fast  Fourie'  transform,  removing  certain  fre¬ 
quencies.  and  then  transforming  back  to  real  space  as  this  results  in  a  periodic  convolution  with 
spurious  end  point  effects.  The  aperiodic  fast  Fourier  transform  is  performed  as  follows.  The 
discrete  Fourier  transform  of  ?,(i)  (denoted  by  ij(*)]: 


NA  1  Zxijk/N.  ... 
r)(k)  =  i  e  1  rp 

j-0 


(A2) 


m,m;  iiw  ppt  aiBori'hni  **  n>-  *  *«■» « *  <c°^.v n*.,  t%5l 

For  a  filtei  witli  2 N  i  1  weights  we  take  N{  >  N  and  define 

C(n  i  N)  n  =  0 . 2N  i  1 

0  n  >  (2/V  +  1) 


C,(n) 


(A3) 


cet^vr;:,i,'discrete;°urier  ,ransf°,m  °f  ~ci(n>  *>y  cm  and  ^  ^  ^m  m)  = 

L  l(k)  W'-  Tll«  discrete  Fourier  transform  of  H(k)  is  given  by 


1  Nl~l 
H(n)  -  L  v 

N\  k=0 


-2irink/N 


m 


(A4) 


and  is  equal  to  the  aperiodic  convolution 


N  -l 

H(n)  =  %{n 

m=0 


m)C(m) 


(A5) 


where  yi)  is  periodic;  i.e.,  r,p(i  ,  qN ,)  =  ^),  ,  an  integer.  Because  of  the  zeros  in  C(j),  how- 
evei,  the  aperiodic  convolution  values  q{(i)  defined  in  eq  2  are  given  by  77.(1)  H(n  f  Af)  for 

!  ~  f  l  Tllus>/Vl  -  2N  filtered  points  are  obtained.  Note  that  if  the  convolution 

length,  2 N  t  1,  is  nearly  equal  to  N  j  then  only  a  few  points  will  be  accurately  filtered. 
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Clearly  the  above  procedure  may  be  used  for  filtering  a  long  record  by  operating  with  the  FFT 
algorithm  on  overlapping  sections  of  record  with  the  overlap  depending  on  the  convolution  length 

In  our  application  the  FFT  aperiodic  technique  (using  N  =  2048)  is  about  five  times  faster  than 
conventional  methods. 


